Model-based Control of Upper Extremity Human-Robot Rehabilitation Systems

Stroke rehabilitation technologies have focused on reducing treatment cost while improving effectiveness. Rehabilitation robots are generally developed for home and clinical usage to: 1) deliver repetitive and stimulating practice to post-stroke patients, 2) minimize therapist interventions, and 3) increase the number of patients per therapist, thereby decreasing the associated cost. The control of rehabilitation robots is often limited to black- or gray-box approaches; thus, safety issues regarding the human-robot interaction are not easily considered. Furthermore, despite numerous studies of control strategies for rehabilitation, there are very few rehabilitation robots in which the tasks are implemented using optimal control theory. Optimal controllers using physics-based models have the potential to overcome these issues. This thesis presents advanced impedance- and model-based controllers for an end-effector-based upper extremity stroke rehabilitation robot. The final goal is to implement a biomechanically-plausible real-time nonlinear model predictive control for the studied rehabilitation system. The real-time term indicates that the controller computations finish within the sampling frequency time. This control structure, along with advanced impedance-based controllers, can be applied to any human-environment interactions. This makes them promising tools for different types of assistive devices, exoskeletons, active prostheses and orthoses, and exercise equipment. In this thesis, a high-fidelity biomechatronic model of the human-robot interaction is developed. The rehabilitation robot is a 2 degree-of-freedom parallelogram linkage with joint friction and backlash, and nonlinear dynamics. The mechatronic model of the robot with relatively accurate identified dynamic parameters is used in the human-robot interaction plant. Different musculoskeletal upper extremity, biomechanic, models are used to model human body motions while interacting with the rehabilitation robot model. Human-robot interaction models are recruited for model-in-loop simulations, thereby tuning the developed controllers in a structured resolution. The interaction models are optimized for real-time simulations. Thus, they are also used within the model-based control structures to provide biofeedback during a rehabilitation therapy. In robotic rehabilitation, because of physical interaction of the patient with a mechanical device, safety is a fundamental element in the design of a controller. Thus, impedance-based assistance is commonly used for robotic rehabilitation. One of our objectives is to achieve a reliable and real-time implementable controller. In our definition, a reliable controller is capable of handling variable exercises and admittance interactions. The controller should reduce therapist intervention and improve the quality of the rehabilitation. Hence, we develop advanced impedance-based assistance controllers for the rehabilitation robot. Overall, two types of impedance-based (i.e., hybrid force-impedance and optimal impedance) controllers are developed and tuned using model-in-loop simulations. Their performances are assessed using simulations and/or experiments. Furthermore, their drawbacks are discussed and possible methods for their improvements are proposed. In contrast to black/gray-box controllers, a physics-based model can leverage the inherent dynamics of the system and facilitate implementation of special control techniques, which can optimize a specific performance criterion while meeting stringent system constraints. Thus, we present model-based controllers for the upper extremity rehabilitation robot using our developed musculoskeletal models. Two types of model-based controllers (i.e., nonlinear model predictive control using external 3-dimensional musculoskeletal model or internal 2-dimensional musculoskeletal model) are proposed. Their performances are evaluated in simulations and/or experiments. The biomechanically-plausible nonlinear model predictive control using internal 2-dimensional musculoskeletal model predicts muscular activities of the human subject and provides optimal assistance in real-time experiments, thereby conforming to our final goal for this project

Predictive Simulation of Reaching Moving Targets Using Nonlinear Model Predictive Control

This Document is Protected by copyright and was first published by Frontiers. All rights reserved. it is reproduced with permission.This article investigates the application of optimal feedback control to trajectory planning in voluntary human arm movements. A nonlinear model predictive controller (NMPC) with a finite prediction horizon was used as the optimal feedback controller to predict the hand trajectory planning and execution of planar reaching tasks. The NMPC is completely predictive, and motion tracking or electromyography data are not required to obtain the limb trajectories. To present this concept, a two degree of freedom musculoskeletal planar arm model actuated by three pairs of antagonist muscles was used to simulate the human arm dynamics. This study is based on the assumption that the nervous system minimizes the muscular effort during goal-directed movements. The effects of prediction horizon length on the trajectory, velocity profile, and muscle activities of a reaching task are presented. The NMPC predictions of the hand trajectory to reach fixed and moving targets are in good agreement with the trajectories found by dynamic optimization and those from experiments. However, the hand velocity and muscle activations predicted by NMPC did not agree as well with experiments or with those found from dynamic optimization.The authors would like to thank the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canada Research Chairs program for financial support of this research

On the Relationship Between Muscle Synergies and Redundant Degrees of Freedom in Musculoskeletal Systems

It has been suggested that the human nervous system controls motions in the task (or operational) space. However, little attention has been given to the separation of the control of the task-related and task-irrelevant degrees of freedom.Aim: We investigate how muscle synergies may be used to separately control the task-related and redundant degrees of freedom in a computational model.Approach: We generalize an existing motor control model, and assume that the task and redundant spaces have orthogonal basis vectors. This assumption originates from observations that the human nervous system tightly controls the task-related variables, and leaves the rest uncontrolled. In other words, controlling the variables in one space does not affect the other space; thus, the actuations must be orthogonal in the two spaces. We implemented this assumption in the model by selecting muscle synergies that produce force vectors with orthogonal directions in the task and redundant spaces.Findings: Our experimental results show that the orthogonality assumption performs well in reconstructing the muscle activities from the measured kinematics/dynamics in the task and redundant spaces. Specifically, we found that approximately 70% of the variation in the measured muscle activity can be captured with the orthogonality assumption, while allowing efficient separation of the control in the two spaces.Implications: The developed motor control model is a viable tool in real-time simulations of musculoskeletal systems, as well as model-based control of bio-mechatronic systems, where a computationally efficient representation of the human motion controller is needed

Configuration-Dependent Optimal Impedance Control of an Upper Extremity Stroke Rehabilitation Manipulandum

Robots are becoming a popular means of rehabilitation since they can decrease the laborious work of a therapist, and associated costs, and provide well-controlled repeatable tasks. Many researchers have postulated that human motor control can be mathematically represented using optimal control theories, whereby some cost function is effectively maximized or minimized. However, such abilities are compromised in stroke patients. In this study, to promote rehabilitation of the stroke patient, a rehabilitation robot has been developed using optimal control theory. Despite numerous studies of control strategies for rehabilitation, there is a limited number of rehabilitation robots using optimal control theory. The main idea of this work is to show that impedance control gains cannot be kept constant for optimal performance of the robot using a feedback linearization approach. Hence, a general method for the real-time and optimal impedance control of an end-effector-based rehabilitation robot is proposed. The controller is developed for a 2 degree-of-freedom upper extremity stroke rehabilitation robot, and compared to a feedback linearization approach that uses the standard optimal impedance derived from covariance propagation equations. The new method will assign optimal impedance gains at each configuration of the robot while performing a rehabilitation task. The proposed controller is a linear quadratic regulator mapped from the operational space to the joint space. Parameters of the two controllers have been tuned using a unified biomechatronic model of the human and robot. The performances of the controllers were compared while operating the robot under four conditions of human movements (impaired, healthy, delayed, and time-advanced) along a reference trajectory, both in simulations and experiments. Despite the idealized and approximate nature of the human-robot model, the proposed controller worked well in experiments. Simulation and experimental results with the two controllers showed that, compared to the standard optimal controller, the rehabilitation system with the proposed optimal controller is assisting more in the active-assist therapy while resisting in active-constrained case. Furthermore, in passive therapy, the proposed optimal controller maintains the position error and interaction forces in safer regions. This is the result of updating the impedance in the operational space using a linear time-variant impedance model